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Abstract 

A conceptual  design  methodology  to  predict  the  wave  drag  of  a transonic  wing  for  use  within 
multisdisciplinary  aircraft  design  was  developed.  To  achieve  this,  a database  of  cross  section 
designs  optimized  with  respect  to  total  drag  was  set  up  varying  the  design  parameters  Ma , t/c,  cl 
and  Re.  Mathematical  formulations  for  the  aerodynamic  cross  section  characteristics  total  drag, 
viscous  drag  and  the  local  shock  location  were  derived  from  the  database  as  functions  of  the  design 
parameters.  The  cross  section  wave  drag  was  then  derived  using  these  formulations.  A locally 
infinite  swept  wing  is  assumed  and  simple  sweep  theory  using  the  shock  sweep  angle  is  used  to 
transform  the  wave  drag.  The  wave  drag  of  a 3-D  wing  is  predicted  summing  locally  infinite  swept 
wing  sections  in  spanwise  direction.  The  achieved  drag  prediciton  is  accurate  enough  for  use  within 
conceptual  aircraft  design  and  predicts  well  the  trends  in  wave  drag  development  as  a function  of 
the  design  parameters  Ma:  t/c:  cl , Re  and  the  wing  planform. 


Introduction 

Within  the  scope  of  multidisciplinary  conceptual 
aircraft  design  [1]  [2],  a transonic  aerodynamic 
module  to  predict  transonic  effects  of  a wing 
is  essential.  Typically,  about  300  parameters 
and  50  design  variables  and  constraints  are  used 
to  describe  the  complete  configuration,  whereas 
about  150  parameters  and  15  design  variables 
of  these  are  used  for  aerodynamic  purposes  [3]. 
The  objective  of  the  presented  methodology  is 
to  describe  transonic  aerodynamic  effects  with 
a limited  number  of  parameters  and  predict  the 
correct  trends  within  conceptual  aircraft  design. 
Processing  time  becomes  critical  for  application 
of  the  methodology  within  multidisciplinary  de- 
sign automation  and  optimization.  For  this 
reason,  time  intensive  codes  (e.g.  3-D  Navier 
Stokes,  3-D  Euler)  are  not  appropriate. 
Therefor  a fast  conceptual  design  methodology 
to  predict  the  wave  drag  of  a transonic  wing  as 
a function  of  the  design  parameters  Mz,  t/c,  c/,, 
Re  and  the  wing  planform  will  be  described. 


Methodology 

Transonic  aerodynamic  characteristics  of  a wing 
are  a result  of  the  wing  planform,  the  spanwise 
distribution  of  lift,  thickness  and  local  flow  data. 
These  parameters  shall  be  used  as  input  data 
for  a method  to  predict  aerodynamic  effects  of 
a transonic  wing. 

As  a first  step,  a database  of  cross  sec- 
tion designs  was  set  up  using  numerical  opti- 
mization routines.  To  achieve  this,  Daimler- 
Chrysler Aerospace  Airbus  (DA)  inhouse  codes 
were  linked  to  the  Synaps  optimization  tool 
PointerPro™  [4]  to  set  up  the  database.  The 
aerodynamic  codes  used  were  VICWA  [5]  (Full 
potential  code  with  finite-difference  boundary 
layer  calculation)  and  XLS6  (Full  potential  code 
with  integral  boundary  layer  calculation).  The 
database  contains  about  240  cross  section  de- 
signs optimized  with  respect  to  total  drag,  called 
optimized  performance  designs.  Typical  opti- 
mized performance  cross  section  designs  of  the 
database  are  shown  in  Fig.  1. 
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All  calculations  are  2D,  without  sweep  angle. 
The  design  parameters  to  set  up  the  database 
of  optimized  performance  cross  section  designs 
have  been  varied  within  the  following  ranges: 

0.7  < Ma2d  < 0.8 
0-5  < ci2d  < 0.9 
0.05  < t/c2d  < 0.16 

For  a 3D  wing  the  corresponding  Ma  range 
for  quarter  chord  sweep  angles  of  about  30°  is 
0.8  < Ma^d  < 0.95.  The  Reynolds  number  Re 
was  varied  from  1-106  to  80 -106. 

Using  this  database  of  optimized  performance 
cross  section  designs,  analytical  approximations 
for  the  relationships  between  aerodynamic  cross 
section  characteristics  and  the  design  parame- 
ters Ma,  t/c,  cl  and  Re  were  derived  from  the 
database  using  appropriate  methods.  Formulas 
were  found,  ensuring  physically  plausible  trends 
towards  the  borders  of  design  parameter  ranges. 
The  idea  for  a methodology  to  predict  the  2D 
total  profile  drag  as  a function  of  the  design  pa- 
rameters Ma,  t/c,  cl  at  constant  Re  originates 
from  A.  Van  der  Velden. 

A refined  methodology  to  predict  the  2D  cross 
section  total  drag  and  a methodology  for  the 
profile  viscous  drag  as  a function  of  the  design 
parameters  Ma,  t/c,  cl  and  Re  were  established. 

Methodologies  needed  for  the  2D  wave  drag: 

• Total  profile  drag 

• Profile  viscous  drag 

To  predict  the  wave  drag  of  a transonic  3D  wing, 
locally  infinite  swept  wing  sections  are  assumed 
which  are  integrated  in  spanwise  direction.  Sim- 
ple sweep  theory  using  the  local  shock  sweep  an- 
gle is  used  to  transform  the  cross  section  wave 
drag.  To  achieve  this,  methodologies  for  the  fol- 
lowing aerodynamic  characteristics  as  functions 
of  the  design  parameters  were  developed. 

Methodologies  needed  for  the  3D  wave  drag: 

• Profile  wave  drag 

• Local  shock  sweep  angle 

The  previously  mentioned  methods  for  the  2D 
total  profile  drag  and  the  profile  viscous  drag 
were  used  to  derive  the  profile  wave  drag.  A 


methodology  for  the  local  shock  sweep  angle  was 
developed  in  previous  work  [6]. 

The  methodologies  are  used  within  multidisci- 
plinary design  automation  to  predict  the  min- 
imum wave  drag,  the  minimum  total  drag  and 
trends  in  drag  development  of  a transonic  wing 
as  a function  of  local  values  of  Ma,  c/,,  t/c,  Re 
and  the  wing  planform. 

The  3D  design  paramter  ranges  for  the  use  of 
the  methodologies  are  as  follows: 

0.7  < Ma  < 0.95 
0.1  < cl  < 1.3 
0.03  < t/c  < 0.18 


Analysis 

The  total  drag  output  from  the  transonic  code 
of  the  optimized  performance  designs  in  the 
database  was  analyzed  to  determine  the  rela- 
tionship between  cross  section  characteristics 
and  design  parameters.  A drag  break  down  into 
wave  drag  due  to  lift  and  viscous  effects  is  as- 
sumed to  be  adequately  correct  for  use  within 
conceptual  design.  Both  wave  drag  due  to  lift 
and  viscous  drag  are  assumed  to  be  functions  of 
the  design  parameters  Ma,  t/c,  cl  and  Reynolds 
number  Re. 


Total  profile  drag 

To  describe  the  total  profile  drag  of  the  opti- 
mized performance  designs  in  the  database,  an 
appropriate  function  was  chosen.  As  a first  es- 
timate a third  order  polynomial  (1)  was  used 
with  mixed  terms  in  the  design  parameters  Ma, 
Cl  and  t/c  combined  with  an  exponential  con- 
tribution of  Ma,  cl  and  t/c  . 
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where 

Xi  = f(Ma ) 

X2  = /(cl) 

^3  = f(t/c). 

For  better  model  handling,  x\ , ^2  and  2:3  are 
linear  functions  of  the  design  parameters  M2, 
c^  and  t/c,  respectively. 

For  a specific  cross  section  n in  the  database, 
the  local  quadratic  deviation  Eiocn  of  the  to- 
tal drag  model  compared  to  the  total  drag  in 
the  database  is  described  in  equation  (2).  Here 
{cD)totdbn  and  ( cD)totmodn  are  the  total  drag  of 
the  cross  section  n in  the  database  and  of  the 
model  following  equation  (1)  respectively. 

-&locn  = (|  {cD)totdbn  — {pD)totmodn  |)  (2) 


database  total  drag  is  shown  as  dashed  lines  and 
the  model  total  drag  as  solid  lines.  Similar  plots 
at  cl  — 0.7  and  cl  = 0.9  are  shown  in  Fig.  3 
and  Fig.  4 respectively. 

The  finally  derived  formula  of  the  total  drag 
model  principally  consists  of  an  absolute  term, 
linear  terms  in  cl  and  t/c  and  functions  of 
(t/c)2.  Third  order  terms  in  the  polynomial 
contribution  were  not  relevant.  An  exponential 
function  with  linear  terms  in  Ma  and  cl  and  a 
quadratic  term  in  t/c  was  used  to  describe  the 
total  cross  section  drag  for  a specific  Reynolds 
number. 

Ad  tot  = Ptot(cL,  t/c,  {t/c)2)  + ■■• 

■ ■ • + Etot{Ma,  cL,  {t/c)2)  (4) 

With  Ptot  polynomial 

Etot  exponential  function 


The  mean  global  deviation  Eg ^ of  the  mathe- 
matical model  for  the  total  drag  for  all  N cross 
sections  in  the  database  may  then  be  derived 
using  equation  (3). 
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To  achieve  the  best  fit  possible  for  the  model 
with  the  database,  the  coefficients  and  Bijk 
were  varied  until  the  minimum  of  the  global  de- 
viation Egiob  was  reached.  Due  to  the  resulting 
nonlinear  equation  system,  a numeric  optimizer 
was  used  to  find  the  coefficients  A ijk  and  Bijk 
in  equation  (1). 

Primarily  a physically  correct  reproduction  of 
trends  in  total  drag  development  was  intended. 
The  mean  deviation  of  the  model  compared 
to  the  total  drag  of  the  cross  sections  in  the 
database  is  less  than  7.5%.  An  exact  agreement 
of  the  total  drag  determined  from  the  database 
with  the  total  drag  calculated  using  equation  (1) 
was  not  required  for  multidisciplinary  design. 
The  total  drag  of  the  optimized  cross  sections 
in  the  database  compared  to  the  total  drag 
model  for  constant  Reynolds  number  at  cl  = 0.5 
is  shown  in  Fig.  2.  The  total  drag  coefficient 
(cD)tot  is  plotted  as  a function  of  t/c.  The 
Mach  number  is  a curve  parameter  and  is  in- 
creasing from  bottom  right  to  top  left.  The 


The  total  drag  model  for  constant  Reynolds 
number  at  cl  = 0.5  is  shown  in  Fig.  5.  The  total 
drag  coefficient  ( co)tot  is  plotted  as  a function 
of  t/c  and  Mach  number  Ma.  Mach  number  is 
increasing  from  front  to  back  and  thickness  is 
increasing  from  left  to  right. 

The  total  drag  model  predicts  the  correct  trends 
in  total  drag  development  as  a function  of  the 
design  parameters  Ma,  t/c,  cl  and  Reynolds 
number  Re.  The  behaviour  near  the  borders  of 
the  parameter  ranges  is  physically  plausible. 

In  order  to  extract  the  profile  wave  drag  from 
the  profile  total  drag,  the  profile  viscous  drag  as 
a function  of  the  design  parameters  was  needed. 


Profile  viscous  drag 

As  for  the  total  drag  model  a similar  approach 
for  the  viscous  drag  model  was  used  to  find  an 
analytical  relationship  of  the  viscous  drag  de- 
pendent on  the  design  parameters  t/c,  cl  and 
Re.  The  optimized  performance  cross  section 
designs  in  the  database  were  used  for  the  calcu- 
lation of  viscous  drag  at  Reynolds  numbers  in 
the  range  from  Re  = 1 • 106  to  Re  = 80  • 106. 

A second  order  polynomial  with  mixed  terms  in 
the  design  parameters  cl  and  t/c  was  used  as  a 
starting  point  for  a specific  Reynolds  number. 
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2 

cDvisc\Re  ~ ^ A?fc  x2  X3  (^) 

j,k=0 
j+h<  2 


Again,  the  coefficients  of  the  mathematical  func- 
tion to  describe  the  viscous  drag  were  found 
using  numerical  optimizers.  The  viscous  drag 
model  describing  the  viscous  drag  of  the  cross 
section  designs  in  the  database  as  a function  of 
the  design  parameters  is  shown  in  equation  (6). 
Note  there  is  no  contribution  of  Mach  number 
to  the  profile  viscous  drag. 


CDvisc  = Pvisc(cL,  t/c,  ( t/c )2)  (6) 

with 

Pvisc  polynomial 

The  coefficients  Ajt~  in  equation  (5)  are  functions 
of  the  Reynolds  number  Re.  Logarithmic  terms 
or  combinations  of  logarithmic  and  exponential 
terms  are  used,  dependent  on  the  coefficient  Aj^. 

A00  = s0  e(lns^) 

Aoi  = go + ln(<?ii2e) 

Aio  = r0  + In  (nife) 

Aq2  = po+PiRe 


The  parameters  pi , q r*  and  Si  in  (7)  were 
calculated  as  previously  described  using  a nu- 
merical optimizer  to  find  a best  possible  fit 
for  the  viscous  drag  derived  from  the  database. 
The  Reynolds  numbers  were  in  the  range  from 
Re  = 1 * 106  to  Re  = 80  * 106. 

The  viscous  drag  of  the  optimized  cross  sections 
in  the  database  compared  to  the  viscous  drag 
model  for  constant  Reynolds  number  at  cl  = 0.5 
is  shown  in  Fig.  6.  The  viscous  drag  coefficient 
( cd)v%sc  is  plotted  as  a function  of  t/c.  The 
Mach  number  is  a curve  parameter  and  there  is 
no  variation  of  viscous  drag  with  Mach  number. 
The  database  viscous  drag  is  shown  as  dashed 
lines  and  the  model  viscous  drag  as  solid  lines. 
A similar  plot  at  ci  — 0.7  is  shown  in  Fig.  7. 
The  viscous  drag  model  for  constant  Reynolds 
number  at  cl  = 0.5  is  shown  in  Fig.  8.  The  vis- 
cous drag  coefficient  (c]j)vi$c  is  plotted  as  a func- 
tion of  t/c  and  Mach  number  Ma.  Mach  number 


is  increasing  from  front  to  back  and  thickness  is 
increasing  from  left  to  right. 

The  viscous  drag  model  predicts  the  correct 
trends  in  viscous  drag  development  as  a function 
of  the  design  parameters  cx,,  t/c  and  Re.  Addi- 
tional drag  due  to  possible  separation  at  the  up- 
per surface  trailing  edge  at  off-design  conditions 
(towards  higher  t/c  at  high  cl)  is  not  described 
with  the  viscous  drag  model  (Fig.  9).  However, 
these  conditions  are  not  relevant  for  transonic 
aircraft  sizing. 


Profile  wave  drag 

Using  the  profile  total  drag  model  and  the  pro- 
file viscous  drag  model  described  in  the  previ- 
ous sections,  a wave  drag  model  can  be  estab- 
lished using  the  assumption  about  drag  break 
down  in  the  section  ’Analysis’.  Subtracting  the 
viscous  drag  model  at  a certain  Reynolds  num- 
ber from  the  total  drag  model  set  up  for  the 
same  Reynolds  number,  the  resulting  model  de- 
scribes the  wave  drag  contribution  of  the  opti- 
mized cross  section  designs  in  the  database  as  a 
function  of  Ma,  t/c,  cl  and  Re.  The  wave  drag 
model  comprises  the  wave  drag  due  to  shock 
waves  and  pressure  drag  due  to  thickening  of 
the  boundary  layer  caused  by  shock  boundary 
layer  interaction. 

Combining  the  wave  drag  model  with  the  vis- 
cous drag  model,  the  total  drag  of  the  cross  sec- 
tions can  be  predicted  as  a function  of  the  design 
parameters  Ma,  t/c,  cl  and  Re. 

The  wave  drag  model  for  constant  Reynolds 
number  at  cl  = 0.5  is  shown  in  Fig.  10.  The 
wave  drag  coefficient  (co)wave  is  plotted  as  a 
function  of  t/c  and  Mach  number  Ma.  Mach 
number  is  increasing  from  front  to  back  and 
thickness  is  increasing  from  left  to  right. 

Wing  wave  drag 

Based  on  the  assumption  of  an  isobaric  wing  de- 
sign concept,  the  wave  drag  can  locally  be  cal- 
culated using  simple  sweep  theory  based  on  the 
local  shock  sweep  angle  (Fig.  11).  The  3D  tran- 
sonic characteristics  of  a wing  may  then  be  pre- 
dicted using  the  2D  models,  derived  from  the 
database.  The  wing  wave  drag  can  be  described 
as  a function  of  the  local  values  of  Ma,  ex,,  t/c, 
Re , the  local  shock  sweep  angle  and  the  wing 

planform. 
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The  local  shock  sweep  angle  is  a function  of  the 
local  shock  location  of  the  upper  surface  shock 
and  the  wing  planform.  An  analytical  model  of 
the  local  shock  location  as  a function  of  the  de- 
sign parameters  was  derived  from  the  database 
simimlar  to  the  drag  models  described  previ- 
ously. The  shock  location  of  the  upper  sur- 
face shock  of  the  optimized  cross  sections  in 
the  database  compared  to  the  shock  location 
model  at  cl  = 0.7  is  shown  in  Fig.  12.  The 
shock  location  ( x/c)s  is  plotted  as  a function  of 
t/c.  The  Mach  number  is  a curve  parameter 
and  is  increasing  from  bottom  right  to  top  left. 
The  database  shock  location  is  shown  as  dashed 
lines  and  the  model  shock  location  as  solid  lines. 
The  shock  location  model  at  cl  = 0.7  is  shown 
in  Fig.  13.  The  shock  location  (x/c)s  is  plot- 
ted as  a function  of  t/c  and  Mach  number  Ma. 
Mach  number  is  increasing  from  front  to  back 
and  thickness  is  increasing  from  left  to  right. 
More  detailed  information  about  the  shock  lo- 
cation model  may  be  found  in  [6]. 

Using  the  viscous  drag  model,  the  wave  drag 
model  and  the  methodology  described  within 
this  paper,  the  minimum  total  drag  of  a 3D 
wing  for  a specific  design  may  be  predicted  and 
a drag  break  down  into  wing  viscous  drag  and 
wing  wave  drag  is  possible. 

The  wing  drag  prediction  solely  is  a function 
of  the  local  design  parameters  Ma,  c^,  t/c,  Re 
and  the  wing  planform.  With  the  drag  mod- 
els being  derived  from  cross  section  designs  op- 
timized with  respect  to  minimum  total  drag, 
the  methodology  predicts  the  minimum  possi- 
ble drag  for  a specific  wing  design.  With  the 
assumption  of  an  isobaric  wing  design  concept, 
a good  prediction  for  the  outboard  wing  and  a 
first  good  estimate  for  the  inboard  wing  may  be 
achieved  using  the  methodology.  However,  the 
more  an  isobaric  wing  design  concept  is  realized 
in  actual  wing  design,  the  better  the  wing  drag 
prediction  which  can  be  achieved  using  the  de- 
scribed methodology. 


drag  development  of  a wing  can  be  predicted  as 
a function  of  the  local  design  parameters  Ma,  cl, 
i/c,  Re  and  the  wing  planform.  As  a result  of  the 
analytical  description  of  the  drag  models  used 
within  the  methodology,  the  processing  time  for 
a wing  drag  prediction  is  only  in  the  order  of 
seconds.  Although  the  methodology  is  based  on 
2D  cross  section  designs  optimized  with  respect 
to  minimum  drag,  the  predicted  minimum  wave 
drag  of  a 3D  wing  compares  well  with  the  wing 
wave  drag  of  current  design  studies.  Future  work 
will  consider  the  implementation  of  effects  of  the 
fuselage  and  engine  installation. 
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Conclusion 

A conceptual  design  methodology  was  described 
to  predict  the  wave  drag,  viscous  drag  and  total 
drag  of  a transonic  wing.  The  methodology  is  in- 
tended for  use  within  multidisciplinary  concep- 
tual design.  The  wave  drag  and  trends  in  wave 
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(a)  Optimized  performance  design, 
Ma  = 0.72,  cL  = 0.7,  (t/c)  = 0.11 


(b)  Optimized  performance  design, 
Ma  = 0.76,  cL  = 0.5,  (t/c)  = 0.10 
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Fig.  2:  Database  total  drag  CDtotdb  versus  formula- 
tion  cDtotmod{Ma,  cL,  t/c)  at  cL  = 0.5 
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Fig.  3:  Database  total  drag  CDtotdb  versus  formula- 
tion  cDtotmod(Ma,  cL,  t/c)  at  c/,  = 0.7 
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Fig.  4:  Database  total  drag  CDtotdb  versus  formula- 
tion  cDtotmod(Ma , cL,  t/c)  at  cL  = 0.9 


Fig.  1:  Typical  cross  section  designs  of  the 
database 
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Fig.  11:  Local  shock  sweep  angle,  isobaric  wing  de- 
sign concept 


Fig.  12:  Database  shock  location  ( x/c)sdb  versus 
mathematical  model  (x/ c)smod(Ma , t/c)  at 
cL  = 0.7 


Fig.  13:  Shock  location  model  (%  / c)  s mod(Ma , t/c) 
at  cl  = 0.7 


